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We study shallow moist Rayleigh-Benard convection in the Boussinesq approximation in three- 
dimensional direct numerical simulations. The thermodynamics of phase changes is approximated 
by a piecewise linear equation of state close to the phase boundary. The impact of phase changes on 
the turbulent fluctuations and the transfer of buoyancy through the layer is discussed as a function 
of the Rayleigh number and the ability to form liquid water. The enhanced buoyancy flux due to 
phase changes is compared with dry convection reference cases and related to the cloud cover in the 
convection layer. This study indicates that the moist Rayleigh-Benard problem offers a practical 
framework for the development and evaluation of parameterizations for atmospheric convection. 
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I. INTRODUCTION 

Moist thermal convection combines turbulent convection with phase changes and latent heat release. It is ubiquitous 
throughout the atmosphere of the Earth (Heintzenberg & Charlson 2009). When a parcel of air rises in convective 
motion, it expands adiabatically. As a consequence, its temperature and pressure drop and at some point during its 
ascent the air parcel becomes saturated. Once water condenses, a cloud is formed. The range of spatial and temporal 
scales in the convective turbulent motion varies widely, from a few hundred meters in isolated cumulus clouds to 
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several thousands of kilometers in midlatitudes storm systems. 
. ■ Despite its enormous importance, the small-scale structure and statistics of moist convective turbulence has been 
studied relatively little compared to its dry convection counterpart. The reason for this gap is that turbulent convection 
j^H *. in moist air includes the complex nonlinear thermodynamics of phase changes in addition to the turbulent motion 
(Stevens 2005; Pauluis 2008). The associated latent heat release provides a rapidly changing local source of buoyant 
i— i, motion, so that moist convection is characterized by a complex interaction between dynamics and thermodynamics. 
One approach to this problem is to express the buoyancy of a parcel of moist air as function of its entropy, pressure 
and total water content. In such framework, phase changes can be treated implicitly, and lead to discontinuities 
ly-j ■ of the partial derivatives in the equation of state at the saturation point (Emanuel 1994). While moist convection 
j remains poorly understood, significant progress has been made in the last decade in understanding the global and 
7-H i local mechanisms of turbulent heat transfer in dry convection (for a comprehensive review see Ahlers et al. 2009). In 
I/"") ' this work, we aim at transferring some of the numerical analysis concepts from the well-investigated dry convection 
^vq case, such as studies of the Rayleigh number dependence of the heat transfer (Verzicco & Camussi 2003), the flow 
I ' properties in the cell (van Reeuwijk et al. 2008) or the small-scale statistics (Emran & Schumacher 2008) to the 
Q\ \ less-explored moist convection case. 

We propose here to take a first step by considering moist convection in the idealized setting of moist Rayleigh- 
| Benard convection with a linearized thermodynamics of phase changes (Pauluis & Schumacher 2010). On the one hand, 
the model is a straightforward extension of numerical studies in dry Rayleigh-Benard convection in the Boussinesq 
approximation (e.g. Schumacher 2009). On the other hand, it is a generalization of a moist convection model which 
was discussed by Bretherton (1987, 1988) for the linear and weakly nonlinear regime and has not been studied ever 
since. Here, we conduct direct numerical simulations of the turbulent nonlinear stage of moist convection. The present 
work reports systematic parameter investigations to understand the effect of phase change on the turbulent transport 
of buoyancy through the shallow layer. We also discuss the dependence of the cloud cover on the physical parameters 
of the model. 

II. MOIST BOUSSINESQ MODEL 

The buoyancy B in atmospheric convection is given by (Emanuel 1994) 
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B(b,q v ,qi,qi,p) = -g , (1) 
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with g being the gravity acceleration, p a mean density, p the pressure, S the entropy and q v , qi, qi the mixing ratios of 
water vapor, liquid water and ice. In the following, we discuss in brief the sequence of simplifications of the equation 
of state that result in a model of shallow non-precipitating moist convection in the Boussinesq approximation - the 
simplest case that goes beyond the well-known dry convection (Pauluis & Schumacher 2010). First, in the Boussinesq 
approximation the pressure variations about a mean hydrostatic profile are omitted when computing the buoyancy 
(Pauluis 2008) and one is left with B(S, q v ,qi,Qi, z). Second, warm clouds are discussed with qi = 0. Third, we assume 
that the air parcels are in local thermodynamic equilibrium, which means that water vapor and condensed water can 
only co-exist at saturation line. This implies that liquid water is formed whenever a relative humidity of 100% is 
exceeded. Furthermore, no rain can fall out in our model. The two remaining mixing ratios are then combined 
to the total water mixing ratio, qr = q v + qi- This assumption also excludes the possibility of supersaturation. 
Condensation in the Earth's atmosphere occurs primarily through heterogenous nucleation caused by a large number 
of cloud condensation nuclei [n ~ 10 9 m~ 3 ). As a consequence, supersaturation rarely exceeds one per cent (Rogers & 
Yau, 1989). In the absence of condensation nuclei in the fluid, homogeneous condensation may result in much larger 
supersaturation, as in a recent experiment by Zhong et al. (2009) where the condensate is formed at the top plate of 
the convection cell. The assumption of local thermodynamic equilibrium has the practical advantage that, once the 
entropy and pressure are known, the total water content can be separated between the vapor and liquid phases. The 
dependencies of the buoyancy are thus reduced to B(S,qr, z). The buoyancy is still a highly nonlinear function of the 
entropy, total water mixing ratio and height. Fourth, we approximate Basa piecewise linear function of the two state 
variables S, qr at each height z around the phase boundary between gas and liquid. The linearization step restricts 
us to a shallow layer since the height variations of thermodynamic quantities have to remain small. It preserves the 
main physical ingredient: the discontinuity of partial derivatives (e.g. the specific heat) at the phase boundary. It also 
allows for an explicit determination of whether an air parcel is saturated or not. Finally, since B is a linear function 
of S and qr, we can introduce two new prognostic buoyancy fields, a dry buoyancy field D (which corresponds to 
a liquid water potential temperature) and a moist buoyancy field M (which corresponds to an equivalent potential 
temperature). They are linear combinations of S and qx- Since the state variables S and qr are adiabatic invariants, 
the two new state variables M and D are also conserved by adiabatic transformations. Consequently, the original 
buoyancy B(S, q v , qi, qi,p) is simplified to B(M, D, z), a linear function of the fields M and D which is given by 

B(x,t) = ma,x(M(x,t),D(x,t) - N 2 z) , (2) 

where N s is the Brunt- Vaisala frequency. This is the saturation condition in our model. 
The dry and moist buoyancy fields can be decomposed in 

D{x, t) = D(z) + D'( X , t)=D + Dh ~ Do z + D '( X , t) (3) 

ti 

M(x, f) = M(z) + M'(x, t) = M + M H- M o z + t) (4) 

The variations about the mean linear profiles of both fields have to vanish at z = and H. Equation ((2} can now be 
transformed into 

B = M(z) + max (M 1 , D' + D(z) - M(z) - N 2 s z) . (5) 

Note that the first term on the right-hand side is horizontally uniform. This implies that it can be balanced by a 
horizontally uniform pressure field given by p(z) — —Mqz ~ [(Mh — Mq) / (2H)]z 2 . We can thus remove the mean 
contribution from the buoyancy field without any loss of generality. A dimensionless version of the equations of 
motion is obtained by defining the characteristic quantities. These are the height of the layer H, the free-fall velocity 
Uf = y/H(Mo — Mh), the time Tf = H/Uf, the characteristic pressure Uj, and the buoyancy difference Mo — Mh- 
The equations together with the decompositions © and @ are given by 



(ill / Pv 

— + (u-V)u = -Vp+ J—V 2 u + B(M,D,z)e z (6) 

V-u = (7) 

^ + (u • V)D' = -j^=^D' + p^u z (8) 
ot y/PrRaM na M 

^ + (u • V)M' = 1 V 2 A/' + u z (9) 
ot \/PrRa M 
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These equations contain three non-dimensional parameters, the Prandtl number Pr, the dry and the moist Rayleigh 
numbers Rar> and Rom 

Pr=», Ra D = H ^- DH \ Ra M = H3 W°- M »l (10) 

K VK VK 

Under most circumstances, the amount of water in the atmosphere decreases with height. This implies that the moist 
Rayleigh number should be larger than the dry Rayleigh number, Rom > Rap. In addition to the three parameters 
explicitly present in equations ©-©j two more parameters are hidden implicitly within the definition (JS|) of the 
buoyancy B which is given in dimensionless form by 



B = max (m', D' + SSD + (\ - jjj^ * - CSAz^j . 



(11) 



The so-called Surface Saturation Deficit SSD and the Condensation in Saturated Ascent CSA are then defined as 

SSD= g° — g° and CSA = . (12) 

M - M H M -M H 

These two new non-dimensional parameters respectively measure how close the lower boundary is to saturation, and 
how much water can condense within the atmospheric layer during an adiabatic ascent of a saturated air parcel. The 
larger CSA, the easier is the formation of liquid water and thus of clouds. When Do — Mq is positive, the air at the 
lower boundary is unsaturated, and Do — Mq is proportional to the "water deficit", i.e. the amount of water vapor 
that must be added to the air parcel to become saturated. A positive Surface Saturation Deficit SSD would occur 
over the continents. For convection over the ocean, the lower boundary is neither saturated nor unsaturated, i.e. 
SSD = 0. It is clear that we can consider a subspace of the five-dimensional parameter space only which is spanned 
in general by Rao, Rclm , Pr, SSD and CSA. Therefore, this study is restricted to Pr = 0.7 and SSD = 0. The 
variation of SSD while keeping the other parameters fixed was discussed already in Pauluis & Schumacher (2010). 

The equations of motion are solved by a pseudospectral scheme with volumetric fast Fourier transformations and 
2/3 de-aliasing in a Cartesian slab with side lengths TH x TH x H . Here T is the aspect ratio of the slab. In lateral 
directions x and y, we apply periodic boundary conditions. In the vertical z direction, we apply free-slip boundary 
conditions, 

u z =D'=M' = and ^ ^ ^ ^ . (13) 

az az 

The boundary conditions, which have also been used in Bretherton (1987, 1988), approximate a situation over an ocean 
surface at the bottom and a temperature inversion at the top. Time-stepping is done by a second-order Runge-Kutta 
scheme. Since both buoyancy fields are linearly unstable, the requirements on mesh resolution and time stepping are 
the same as in dry convection. The additional scalar field and the update of the B increases computational costs by 
20%. Table 1 summarizes the grid resolutions and dimensionless parameter sets which are taken in the direct numerical 
simulations. The spectral resolution does not go below k max r]K = 2.45 for all DNS, where k max is the maximum 
resolved wavenumber and t\k the Kolmogorov length. Technically, we use B' in the momentum equation instead of 
B since the mean contribution is B(z) which can be added to the kinematic pressure, i.e. d z p + B = d z p + B'. For the 
moist runs we distinguish two classes for initial equilibrium configurations - a fully saturated slab which corresponds 
with M(z) > D(z) - N^z (large CSA) and a fully unsaturated slab with M(z) < D(z) - N^z (small CSA). 



III. RESULTS 



A. Buoyancy and velocity fluctuations 



Initially the equilibrium configuration is perturbed infinitesimally and after T/Tf ~ 10 2 the flow is relaxed into 
a fully developed and statistically stationary turbulent state. This is when the statistical analysis is started. As 
indicated in the table, we restrict dependencies of moist convection on the two Rayleigh numbers and parameter 
CSA. Figure [1] shows the mean vertical profile of the buoyancy, (B(z)), which is obtained by taking averages in the z 
planes and over an ensemble of statistically independent snapshots. It is observed, that this profile becomes strongly 
asymmetric for the cases which started with an initially fully saturated equilibrium. These will be the cases where 
phase changes affect the turbulence properties most strongly. 

Table 1 lists the root-mean-square (rms) values of M and u as obtained in the statistically stationary regime. 
Since both buoyancy fields follow a linear advection-diffusion equation, the ratio of the rms fluctuations to the outer 
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Run N x x N y x N z Ra M Ra D CSA Uf Tj T/Tj 



r 


512 x 512 x 65 


9.5 


X 


10 5 


7.0 


X 


10 5 





53 


3 


06 


1 


03 


141 





356 





434 


2* 


512 x 512 x 65 


1.1 


X 


10 6 


7.0 


X 


10 5 





44 


3 


35 





94 


154 





334 





436 


3 


512 x 512 x 65 


1.4 


X 


10 6 


7.0 


X 


10 5 





35 


3 


75 





84 


173 





272 





432 


4 


512 x 512 x 65 


1.9 


X 


10 6 


7.0 


X 


10 6 





26 


4 


36 





72 


368 





225 





431 


5 


512 x 512 x 65 


2.9 


X 


10 6 


7.0 


X 


10 5 





17 


5 


35 





59 


449 





184 





433 


6 


512 x 512 x 65 




- 




7.0 


X 


10 5 





00 


2 


63 


1 


19 


122 





362 




— 


7* 


1024 x 1024 x 129 


9.5 


X 


10" 


7.0 


X 


10 6 





53 


3 


06 


1 


03 


342 





308 





436 


8* 


1024 x 1024 x 129 


1.1 


X 


10 7 


7.0 


X 


10 6 





44 


3 


35 





94 


407 





287 





437 


9 


1024 x 1024 x 129 


1.4 


X 


10 7 


7.0 


X 


10 6 





35 


3 


75 





84 


412 





240 





436 


10 


1024 x 1024 x 129 


1.9 


X 


10 7 


7.0 


X 


10 6 





26 


4 


36 





72 


431 





194 





434 


11 


1024 x 1024 x 129 


2.9 


X 


10 7 


7.0 


X 


10 6 





17 


5 


35 





59 


587 





158 





436 


12 


1024 x 1024 x 129 








7.0 


X 


10 6 





00 


2 


63 


1 


19 


321 





320 






13* 


2048 x 2048 x 257 


1.1 


X 


10 8 


7.0 


X 


10 7 





44 


3 


35 





94 


125 





258 





439 


14 


2048 x 2048 x 257 


1.9 


X 


10 8 


7.0 


X 


10 7 





26 


4 


36 





72 


150 





176 





438 



TABLE I: Parameters of simulation runs: grid resolution, RaM, Ran and CSA. For all runs, Pr — 0.7, r = 8 and SSD — 0. We 
also display the characteristic velocity Uf = y (Mo — Mh)H, the characteristic time scale Tf = H/Uf and the total integration 
time T/Tf. For dry runs 6 and 12, Uf = yj (Do — Dh)H. Furthermore, Urms/Uf with u rms — y/ (u| +u y + u1) x ^ y , z ^ and 
M rms /(M - M H ) with M rms = y/{M' 2 ) x , y ,z,t are shown. Runs that are labeled with an asterisk start out of a completely 
saturated equilibrium, M(z) > D(z) — N^z. 
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FIG. 1: Mean total buoyancy profiles (B(z)) for Runs 1 to 5. All runs have the same amplitude of Mo (which equals Do)- It 
is also indicated for which runs the initial equilibrium solution is completely unsaturated or saturated. Profiles for the series 
with RaD — 7 x 10 6 look qualitatively similar, except that the boundary layer thickness decreased. 



buoyancy difference Mq — Mh should be constant. The velocity rms fluctuations decrease with decreasing CSA. At 
fixed CSA, the fluctuations decrease also with increasing Rayleigh number Rap. This result is also observed in dry 
convection (Verzicco & Camussi 2003). In Fig. [3]we refine this analysis and study the vertical profile of the rms of the 
vertical velocity component, it z , rms (z) = yj (u^(z)) XiVi t. A measure of asymmetry of the profile with respect to the 

midplane z = H/2 can be based on the moments m n = (z — -y)" u zjms (z) dz. If the skewness ms/m^ 2 is larger 
than zero then the vertical velocity fluctuations are enhanced in the upper half of the slab. Figure [2] (a) shows that 
the profile is symmetric for the dry reference run and those with smaller amount of water which can be condensed. 
Asymmetry is observed for CSA > 0.35 which peaks at CSA — 0.35 and decreases again for larger CSA. We will 
show at the end of subsection 3.2 that the asymmetry in the vertical velocity fluctuations is directly coupled to the 
vertical fraction Az/H of the convection layer that is partially saturated and unsaturated. This fraction turns out 
to be largest at CSA = 0.35 for all Rap (see Fig. d](b)). It is also found that isotropy in the velocity fluctuations 
is established to a better degree with increasing Rcld- We conclude that phase changes cause the asymmetry of the 
vertical velocity fluctuations. However, with increasing Rayleigh number and thus Reynolds number the small-scale 
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FIG. 2: Asymmetry of the velocity fluctuations as a function of CSA and Rcld- (a) Skewness ms/n^ of the vertical profiles 
of the root mean square of the vertical velocity component u z . (b) Vertical fraction Az/H of the layer with 0.05 < F(z) < 0.95 
where F(z) is given by (|18|) . 




FIG. 3: Convective part of the buoyancy flux as a function of Rcld, Rdu an d CSA. (a) Relation between the fluxes (u z D(z)), 
(u z M(z)) and {u z B(z)}. Data are for CSA = 0.44 and Ra D = 7.0 x 10 6 . (b) Same as (a) for CSA = 0.35. (c) (u z B{z)) 
normalized by the corresponding dry diffusive buoyancy flux which is given as k(Dq — Dh)/H. Data are for runs 7 to 12. 



turbulence is found at increasingly isotropic conditions which can compensate this trend in parts. 

Of central importance in dry convection is the one-point-correlation between buoyancy (or temperature) and vertical 
velocity, (u z B(z)) (which is equal to (u z B' (z))). It enters the definition of the dimensionless measure of buoyancy 
flux through the layer, the Nusselt number Nu. In the present model, we can define two Nusselt numbers for both 
fields in a standard way, such as Nud(z) — [( u zD(z)) — Kd z (D(z))]/\n(DQ — Dh)/H] for D which is constant and thus 
simply denoted by Nud- Since Nun and Num are normalized with respect to their diffusive fluxes, Nun = Num 
follows which was verified in the simulations. In order to quantify the additional amount of buoyancy transfer, we 
will relate the correlations (u z B(z)) to the dry field in the following. Note that the buoyancy flux (u z B(z)) is tied to 
the correlations (u z D(z)) and (u z M(z)). Moreover, because the partial derivative of the buoyancy with respect to M 
and D are bound by and 1 (see ©), we have automatically that at each level z 

(u z D(z)} < (u z B(z)) < (u z M{z)} . (14) 

The lower bound ((u z D(z)) — (u z B(z))) occurs when a layer is fully unsaturated, while the upper bound ((u z B(z)) = 
(u z M(z))) is achieved in fully saturated layer. This is demonstrated in Figs. |3] (a) and (b). In Fig. [3] (c), we 
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FIG. 4: Vertical profiles of Nud(z) and Nub(z) as given by Eq. (|15[) . Data are for Runs 7 and 12. The additional amount 
of buoyancy, Q is the area between moist and dry reference profiles. The inset shows Q as a function of CSA and Rao (Runs 
1-5 and 7-11). 



normalize the correlation by the dry diffusive buoyancy flux (which would correspond with Nu£> = 1)- It is given 
by k(Dq — Dh)/H. Again, we observe an enhancement of the correlations for the three largest values of CSA. The 
profiles for the two smaller values of CSA collapse almost perfectly with the corresponding dry reference cases in 
both series of DNS. Note that a normalization by the moist diffusive buoyancy flux k(Mq — Mh)/H would result 
in systematic growth of the correlation since an increase of Rcim is in line with a decrease of CSA. Finally, the 
correlations increase as well when the Rayleigh numbers Rao and Ra,M are enhanced at given CSA. 

On the basis of the correlations between buoyancy and vertical velocity and the mean vertical profiles, the additional 
buoyancy flux due to phase changes can be determined. We define a Nusselt number based on the dry diffusive 
buoyancy flux: 

K ' k(D -D h )/H ' V ; 

which is not necessarily constant with height, as can be seen in Fig. [4] Similar behaviour was found by Oresta et al. 
(2009) in bubbly convection with phase changes. The additional buoyancy flux due to phase changes and latent heat 
release can be quantified in terms of the parameter: 

1 ^ 



Q=- Nu B {z)dz~Nu D . (16) 
H Jo 

The upper and lower bounds on the buoyancy flux (fUtj) implies that the enhancement factor Q is itself bound by 
1 < Q < Rclm / R&D- Figure [4] shows Q as a function of CSA and Rao- The sensitivity of Q to CSA at a given value 
of Rao is complex. On the one hand, high value of CSA implies more water and a deeper saturated layer. On the 
other hand, in our experimental set-up with SSD = and constant Rajj, any increase of CSA is connected with a 
decrease of Ram- We observe here a maximum of Q at CSA = 0.44. This is the case where in Fig. [3] the largest 
amplitudes for (u z B(z)) are observed (see dashed line in panel (c)). 



B. Cloud cover 



The phase changes in the convective turbulence are associated with the appearance and disappearance of clouds. 
They are defined as those sites where the liquid water mixing ratio qi(x,t) > 0. Translated into our framework this 
corresponds with 

ffi(x,t) = M(x,t) - [X>(x,t) -N^z] > 0. (17) 

The cloud boundary is given by qi = M — D + N^z = 0. Depending on CSA and both Rayleigh numbers this is a 
simply connected isosurface or a collection of disconnected isosurfaces. The latter case is illustrated in Fig. [S] The 



FIG. 5: Snapshot of the instantaneous cloud distribution (white isosurface) in combination with the updrafts (red isosurface). 
The white isosurface is the cloud boundary with qi = 0. The red isosurface is for u z > 0.23(7/ or 1.5it rms . Data are for Run 
14 with CSA = 0.26. 



white isosurfaces qi = display isolated clouds. They are correlated with strong updrafts as illustrated by the red 
isosurfaces for u z > 0.23C//. Warm air rises up and expands adiabatically such that the temperature decreases and 
condensation sets in. 

Figures|5](a) and (b) display the probability to find clouds at height z in the slab as a function of Ran and CSA in a 
semi-logarithmic plot. The formation of clouds is less probable when Ra^ is increased. Reasons could be the stronger 
filamcntation of the turbulent patches and the decreased velocity fluctuations which are in line with an increase of 
the Reynolds number of the turbulent flow. For CSA = 0.53 and 0.44 the cloud layer is closed for all Rao which is in 
line with P(z\qi > 0) = 1. For CSA = 0.35, 0.26 and 0.17, a broken cloud layer with isolated clouds can be observed. 
While the former could correspond with a stratocumulus-like convection regime, the latter could correspond with a 
cumulus-like regime. Note also that layer remains basically dry for the smallest CSA. 

The presence of clouds is also related to the enhancement of the buoyancy flux. We define a saturation fraction 
F(z) for the buoyancy flux which follows from (|14[) and is given by 

_ (u z B(z)) - (u z D(z)} . . 

[ ' (u z M(z)) - (u z D(z)) ■ {LS> 

This saturation fraction is such that in a fully saturated environment F(z) — 1, while in an unsaturated environment, 
we have F(z) — 0. Figures [5] (c) and (d) replot the same data sets for F(z) and in panel (d) a direct comparison 
with P(z\qi > 0) is provided. On the basis of the data, we can conclude that both measures collapse quite well 
for fully saturated or unsaturated layer, but the case CSA = 0.35 indicates there are some significant departures in 
partially saturated layers. Finally, one can define now the vertical fraction of the layer that is partially saturated 
and unsaturated as 0.05 < F(z) < 0.95 (see Fig. |6] (d)). This fraction is biggest for the runs at CSA — 0.35 as 
shown already in Fig. [3](b). In this case, the asymmetry between saturated ascents and unsaturated descents can 
be expected to be largest. Consequently, the largest asymmetries of the vertical velocity fluctuations can be built 
up. The peak at exactly the same CSA value in both panels of Fig. [3] supports our conclusion and closes the loop 
between buoyancy transfer, cloud cover and vertical flow asymmetry. 



IV. SUMMARY AND CONCLUSIONS 



We have presented a shallow moist convection model with a linear equation of state for the thermodynamics of 
phase changes. This model which contains five dimensionless parameters is discussed in a three-dimensional subspace 
due to fixed Prandtl number and Surface Saturation Deficit (SSD). The most important simplification which reduces 
the complexity is the assumption of a local thermodynamic equilibrium. Several key physical processes, such as the 
formation of precipitation and the existence of supersaturation are thus omitted. The model nevertheless captures the 
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FIG. 6: Measures of the cloud cover, (a) Probability P to find clouds in a layer at height z as a function of Rao- Data are for 
CSA = 0.44 . (b) Same as in (a) for CSA = 0.26. The corresponding dry Rayleigh numbers for both figures are given in the 
legend of the left figure, (c) Function F(z) as a function of Rao and CSA. Line styles correspond with (a), (d) Comparison 
of P and F(z) for runs 1 to 5. Data for run 5 at CSA = 0.17 coincide with the axis. 



fundamental interactions between phase transition and dynamics. Phase changes cause an asymmetry of the vertical 
velocity fluctuations when the amount of water that can be condensed (parameter CSA) is sufficiently large. Similar 
to the dry convection case, the correlations between vertical velocity and buoyancy are used to quantify the amount 
of additional buoyancy flux due to condensation and related latent heat release. Furthermore, this correlation can be 
connected with the cloud cover in the layer. 

The studies in this simplified setting provide thus a basis for possible parameterizations of cloud impact in large- 
scale models. In particular, determining the factors that control cloud fraction is a central issue in climate modeling, 
as small changes in cloud cover can dramatically affect the amount of energy received and emitted by the atmosphere. 
We found here that the Rayleigh number has a direct impact on the cloud cover, which should be a cause of concern, as 
the Rayleigh number in our numerical simulations (Rao = 7 x 10 7 ) is significantly smaller than its typical atmospheric 
value, Ra w 10 18 -10 22 . Nevertheless, the idealized moist Rayleigh-Benard convection provides an important test for 
our understanding of clouds and of their sensitivity to environmental parameters. 
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